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INTRODUCTION 

The fascinating physics of the perovskite manganitesi is governed by electrons which 
hop among or are localized on Mn 3+ and Mn 4+ ions. It was pointed out already in 
1951 by Zenera that the strong Hund coupling between e g and t2 g electrons is essential 
in understanding the ferromagnetismi caused by substitution of triply valent La with 
doubly valent elements such as Ca, Sr or Ba in LaMn03. The two- fold degeneracy of e g 
orbitals causes other important effects in the physics of manganites. The existence of 
the orbital degrees of freedom leads to orbital ordering in the insulating phases B Various 
charge and orbital ordered states are experimentally observed]! Another important effect 
is the Jahn- Teller distortion, which causes strong correlations between lattice distortions 
and orbital and magnetic interactions .Em All these effects are entangled and give rise to 
the rich phase diagrams and peculiar transport phenomena of the manganites. Though 
it might be necessary to take into account all of these effects at the same time in order to 
explain experimental results quantitatively, it is useful to extract one of these important 
factors and study its effect in detail, in order to understand the physics of manganites 
in a profound way. In this paper we aim to clarify the effects of the Hund coupling 
by employing simplified models. First we study the ground state phase diagram of a 
doubly degenerate Hubbard model. Our main concern is the effectiveness of the Hund 
coupling on the ferromagnetism. We compare the result in one dimension with that 
in infinite dimensions, and examine common features and differences between them. 



lr To appear in proceedings of the workshop "Physics of Manganites", Michigan State University, July 
26-29, 1998, eds. by T.A. Kaplan and S.D. Mahanti, Plenum Publishing Corporation. 
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The second model we study is a simple ferromagnetic Kondo lattice model or double 
exchange model. We investigate one-particle states in this model using a single-site 
approximation and calculate the electrical resistivity. We treat the localized spins as 
quantum mechanical ones and study the quantum effects on the electronic states. 



Doubly Degenerate Hubbard Model 

We consider in this section a doubly degenerate Hubbard model described by the 
Hamiltonian 



H = -t (4ma c jrna + hx.) + UYn im ^n im i + U' Y nii a n i2a > 

m=l,2 /jjNgN.N. i,m i,cr,cr' 

ct=T,J. 
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where Cj mo - (c] mcr ) denotes the annihilation (creation) operator of the electron at site i 
with orbital m(=l or 2) and spin a. The number operators are denoted by ni ma . Hop- 
pings of electrons are assumed to occur between the same orbitals of nearest neighbor 
sites. In real systems, there are off-diagonal hoppings and also hopping integrals are 
anisotropic for e g orbitals, that is, they are dependent on the directions of hoppings. 
This anisotropy may have an important effect on the orbital and antiferromagnetic or- 
dering in manganitesM We take here, however, the simplest model which can take 
account of the effects of orbital degeneracy and Hund coupling. The interaction terms 
in eq. (|]) originate from the Coulomb interaction between electrons at the same site. 
The last term, which transfers two electrons on one orbital to the other, is often ne- 
glected. But this term should be properly considered, since it enhances local quantum 
fluctuations and the coefficient J' is equal to that of the Hund coupling J if we assume 
the orbital wave functions are real. The interaction parameters satisfy the relation 

U = U' + 2 J (2) 

for e g orbitals, and we assume this relation in the following. 

The Hamiltonian (jl]) produces various magnetic correlations and effective ferromag- 
netic interactions between electrons on different sites. Let us first consider two electrons 
on one site. The spin-triplet states, where two electrons occupy different orbitals, are 
stabilized by the Hund coupling and have the lowest energy U' — J. There are three 
spin-singlet states; one with energy U' + J where two electrons occupy different orbitals, 
and the other two with energies U — J'(= U' + J) and U + J'(= U' + 3J) where elec- 
trons occupy the same orbital. An effective spin interaction between neighboring sites 
is derived from this one-site spectrum in the strong correlation regime (U' > J > t). 
Let us consider two nearest neighbor sites each of which is occupied by a single electron 
in this regime. When two electrons with parallel spins sit on different orbitals, virtual 
hoppings of the electrons between two sites lower the energy by —2t 2 /(U' — J). When 
two electrons have antiparallel spins, the energy is lowered by —2t 2 U/(U 2 — J' 2 ) or 
—2t 2 U'/(U' 2 — J 2 ) depending on whether they are on the same or different orbitals. 
Hence there is an effective interaction between neighboring sites which favors ferro- 
magnetic spin alignment but alternating alignment of orbital degrees of freedom.0 In 
a system with quarter-filled bands, i.e. for n = N e /N = 1, both ferromagnetic long- 
range order (LRO) and alternating orbital order are expected to coexist in the ground 
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state. Here N e and N denote the total number of electrons and sites, respectively. 
The second-order perturbation from the atomic limit (t = 0) leads to the following 
effective Hamiltonian for spin operators Sj and pseudo-spin operators 
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In strongly correlated systems with the filling 1 < n < 2, each lattice site is either 
singly- or doubly-occupied, and electrons hop from doubly-occupied sites to singly- 
occupied ones. Doubly-occupied sites are almost necessarily in spin-triplet states due 
to the Hund coupling and the hopping probability is largest between pairs of sites 
with parallel spins. As a result the kinetic energy is lowered by ferromagnetic spin 
correlations. This mechanism favoring ferromagnetism is quite similar to that in the 
double exchange model of electrons, where electrons interacting with localized spins 
have lower kinetic energy when spins are aligned parallel.^' M In the following we call 
this mechanism which favors ferromagnetism the "double exchange mechanism" even 
when we are not treating localized spins. In the case with less-than-quarter filling 
(n < 1), the "double exchange mechanism" may not work for U' — J 3> t. Nevertheless 
the Hund coupling may lead to ferromagnetism even for n < 1, if t/(U' — J) is not 
too small. The effective ferromagnetic interaction described by eq. (||) between nearest 
neighbor electrons may have a sizable effect and cause metallic ferromagnetism. V; 
Vleck argued that this mechanism may be operative in realizing ferromagnetism in Ni 

Though the mechanism favoring ferromagnetism can be understood qualitatively 
as above, it is far from trivial whether ferromagnetic long-range order occurs in bulk 
systems. In the following, we present a numerical study of the model in one and infinite 
dimensions. 



One-Dimensional Model 



There are rigorous proofs for the ferromagnetic ground state of the one dimensional 
model in strong coupling limits ]ll~§2l These proofs are valid in different limits of strong 
coupling. For the strong Hund coupling case (J — » 00 and U — * 00), existence of 
ferromagnetism is proved for arbitrary U'(> 0) in 1 < n < 2,0 @ and also for < n < 1 
in the special limit J = U' — > 00 and U — > 00 .113 Shen obtained a rather general result 
that the ground state is fully spin-polarized for any n between and 2 except for 1 if 
U = 00, and U'(> 0) and J = J'(> 0) are finite£3 (We note that this result cannot be 
applied naively to our case which assumes the relation (|2|).) 

So far several numerical studies were done, and ferromagnetism was found for den- 
sities near quarter filling.0' &§J These studies were done by diagonalizing relatively 
small systems with sizes up to 12, and size dependence was not studied yet. We thus 
need to study systems with larger sizes and examine size effects to obtain conclusive 
results. We note also that most previous studies did not take into account the J'-term, 
and assumed relations between U, U' and J which differ from eq. 
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We report in the following a study of finite-size chains with up to 16 sites applying 
the exact diagonalization method as well as the density matrix renormalization group 
(DMRG) method. We employ open boundary conditions, since the periodic boundary 
condition causes very large size dependence in one dimension (e.g. even-odd oscilla- 
tions). Remarkably we found little size dependence due to the use of the open boundary 
conditions. Details of this study will be published elsewhereHil 

First, we show the ground-state phase diagram for n = 1 in Fig. la. We obtained 
the ground state with full spin polarization for J ~ U'Qt 5i). The appearance of ferro- 
magnetism in the J < U' region can be well understood with the effective Hamiltonian 
(H). This ferromagnetic ground state has a strong alternating correlation in the or- 
bital degrees of freedom. This is also consistent with the argument from the effective 
Hamiltonian. In the perfect ferromagnetic ground state, the orbital degrees of freedom 
have isotropic (Heisenberg) antiferromagnetic interaction with pseudo-spin and hence 
the alternating correlation decays in a power form. The phase boundary for J < U' ap- 
proaches an asymptote J = aU' with a ~ 0.35 for large U' . The asymptote corresponds 
to the ground-state phase boundary of the effective HamiltonianE3 The paramagnetic 
state for ^ J <LodJ' may have the same properties as the ground state of the SU(4) 
model at J = O.Hc3 The ferromagnetic phase extends to the parameter region J > U', 
as well. Though this region J > U' is unrealistic, it is of interest from the viewpoint of 
triplet superconductivity. An attractive force acts between electrons with parallel spins 
due to the Hund coupling and the present model might have some relevance to exper- 
imental results on organic superconductors.il For J > U' 3> t two electrons are paired 
to form a hard-core boson with spin unity. The perturbation due to the hopping term 
in eq. (1) leads to the effective Hamiltonian for these bosons, which includes hopping, 
repulsion and antiferromagnetic spin interaction between nearest neighbor bosons. The 
effective Hamiltonian does not favor ferromagnetism and indeed the ferromagnetic phase 
does not extend to a region with large J — U' in our numerical calculations. The slope 
of the phase boundary approaches unity for large U'. We note that only the systems 
with N e = even are used in determining the phase diagram. In fact we found a large 
difference in the phase boundaries for J > V between systems with even iV e and odd 
N c . We consider that the results for small odd iV e are strongly affected by the existence 
of an unpaired electron and are not useful for extracting bulk properties. 

Next we show the ground state phase diagram for n = 0.5 in Fig. lb. The ferromag- 
netic phase expands compared to that in the quarter-filled case both for J > U' and 
J < U', and the size dependence is very weak. It is remarkable that ferromagnetism is 
realized for rather weak Hund coupling, that is, J ~ 2t for U' ~ 5. Recently HirschBJ 
argued that the Hund coupling is not effective enough to realize ferromagnetism in sys- 
tems with low density (n < 1) and that ferromagnetic exchange interactions between 
different sites are necessary to explain ferromagnetism in low density systems like Ni. 
From the present results, we expect that a moderate Hund coupling realizes ferromag- 
netism in a bulk system in one dimension. This behavior should be compared with the 
result for infinite dimensions where we could not find ferromagnetism for n < 1. (See 
next section.) 

Finally, as an example of the case with n > 1, we show the result for n = 1.25 in Fig. 
lc. In this case ferromagnetism appears in a wider region than in the quarter filled case 
especially for small J. The lower phase boundary apparently approaches the line J = 
for large U'. This enhanced stability of the ferromagnetic state may be understood as 
the result of the "double exchange mechanism" . On the other hand the phase boundary 
for J > U' is almost same as that for n = 1. 

For all densities we found ferromagnetism on the line J = U' with strong J(^> 5t). 
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Figure 1: Ground-state phase diagrams of the ID doubly degenerate Hubbard model 
for the filling n = 1(a), 0.5(b) and 1.25(c). We assume t to be unity. 



This result is consistent with the rigorous result in the limit J = U' — » oo by Kusakabe 
and AokiJll We note that all the ferromagnetic ground states obtained above are fully 
polarized. Since the ferromagnetic state is fully polarized the spin degrees of freedom 
are completely suppressed. The orbital degrees of freedom in the ground state are 
mapped into the usual spin degrees of freedom in the single-band Hubbard model with 
the interaction parameter U' — J.clJ Then we learn that, for U' — J > 0, the pseudo-spin 
(orbital) correlation function decays with a power law as cos(|i — j\nir) ■ \i— j\~ a . On the 
other hand, for U' — J < 0, it shows an exponential decay but the pair-pair correlation 
function of pseudospin-singlet (spin-triplet) pairs decays with a power law, which is a 
sign of quasi-long-range order of the triplet superconductivity. 

Infinite Dimensional Model 

Next we discuss the model (|l]) on a hypercubic lattice in infinite dimensions .0 We 
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scale the hopping integrals between nearest neighbor sites as t = t/2yd in d-dimensions 
and consider the limit d = oo. Then the density of states (DOS) of each energy band 
has the Gaussian form D(s) = exp(— e 2 /ty/it . We assume t = 1 in the following. In 
this limit we can treat quantum fluctuations completely by taking local interactions into 
account and spatial correlations can be neglected^ The system is described in terms 
of a one-site effective action which is determined self-consistently. Generally one must 
rely on numerical methods to solve the effective action. In this study we approximated 
the action by that of a two-channel impurity model with finite (n s ) number of levels 
in each channel, and solved the impurity model by exact diagonalization. The energy 
levels and mixing parameters of the impurity model were determined self-consistently. 
We searched for ground states which are uniform in space as well as those with two- 
sublattice structures. Numerical calculations were done for n s = 5 or 6. We mostly 
studied the system with n s = 5 and confirmed phase boundaries by using the system 
with n s = 6. We found that the results do not depend much on n s . We studied the 
ground state mainly at the fillings n = 1, 1.2, and 0.8, controlling the chemical potential. 




Figure 2: Ground state phase diagram of the D = oo doubly degenerate Hubbard model 
for n = 1(a) and 1.2(b). 

In the quarter-filling (n = 1) case we found paramagnetic and ferromagnetic ground 
states for < J < U' as is shown in Fig. 2a. Near the phase boundary two solutions 
coexist and their energies cross over. We selected the ground state by comparing ener- 
gies and determined the phase diagram. The paramagnetic state obtained is spatially 
uniform and metallic. On the other hand the ferromagnetic state has a two-sublattice 
structure with alternating orbital order, and is insulating. We found a narrow para- 
magnetic region for J ~ U'. Therefore the ferromagnetic phase seems to be confined 
within the region J < U' , though we did not study the case with J > U' in d = oo. The 
ferromagnetic ground state appears even for J ~ for U' Jo" 6i. At J = the model 
possesses an SU (4) symmetry. We found the coexistence of several ground states in this 
case. 

We show the phase diagram for n = 1.2 in Fig. 2b. At this filling we obtain a metal- 
lic ferromagnetic phase and a metallic paramagnetic one, both of which are spatially 
uniform. Both states have no orbital ordering. The area of the ferromagnetic phase is 
reduced in the phase diagram compared to that of the insulating ferromagnetic phase 
at n = 1. As a hole-doped case, we studied the ground state for n = 0.8. At this filling 
we found only metallic ground states which are uniform in space and could not find any 
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Figure 3: The number density as a function of the chemical potential /i for J7' = 10 and 
J = 4. The flat parts indicate insulating states. 

magnetically ordered phase for U' < 20i. 

We show the variation of the number density as a function of the chemical potential 
for U' = 10 and J = 4 in Fig. ^. For these parameter values we have the paramagnetic 
metallic ground state for < n < 0.82. For n = 1 the ground state is a ferromagnetic 
insulator. The ferromagnetic metal is stable for 1.14 < n < 1.86. The antiferromagnetic 
insulator is realized for n = 2. It is interesting that there are small jumps of n on both 
sides of quarter-filling. One is from n = 0.82 to 1 and the other is n = 1 to 1.14. There 
is another jump close to half-filling, i.e. between n = 1.86 and 2. They imply that phase 
separation occurs for n in these intervals. Occurrence of phase separation was found 
also in the double exchange model£3 

We show the kinetic and interaction energy per site as a function of U' on the line 
J = 0AU' for n = 1.2 in Fig. |j. For these parameters the ground state is ferromagnetic 
for U' 6. The kinetic energy increases linearly with U' for small U' where the ground 
state is paramagnetic. Then it starts to saturate and stays almost constant in the 
ferromagnetic phase. The interaction energy increases linearly with U' but its slope 
decreases as an effect of (local) correlations. The slope slightly increases again in the 
ferromagnetic region and the potential energy is nearly (n — 1)(( 7 7 / — J) = 0.12U' for 
large U' . The above result clearly shows that the ferromagnetism is caused by reduction 
of kinetic energy rather than interaction energy. That means that the "double exchange 
mechanism" is the cause of ferromagnetism for this density. 

We have seen above that the Hund coupling is effective both in one and infinite 
dimensions. Especially for n > 1 the ferromagnetic ground state is realized in quite 
a large parameter region in both dimensions. This result suggests that the "double 
exchange mechanism" is quite effective in realizing the ferromagnetic state for 1 < n < 2. 
We may expect that the situation is similar in two and three dimensions. For n = 1, 
ferromagnetism accompanied by alternating orbital order is realized for J < U'. We 
found that the state is destabilized for weak J in one dimension. This result may 
be understood as lower dimensionality stabilizing the paramagnetic liquid state and 
destabilizing alternating orbital order. (In fact we have only quasi-long-range orbital 
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Figure 4: The kinetic and the potential energy per site vs U' for n = 1.2 on the line 
J = OAU'. 

order in one dimension.) The paramagnetic state for J = is known to be an SU (4) 
singlet state&H 

We found that the phase diagram is strongly dependent on the dimensionality for 
n < 1. Though ferromagnetism is realized in a large parameter region in one dimension, 
we could not find it in infinite dimensions. Although our study does not exhaust the 
whole parameter region, it seems likely that there is no ferromagnetism for n < 1 in 
d = oo. There may be a general tendency for ferromagnetism in a low density system 
to be stabilized in one dimension. In a ID Hubbard model with nearly flat bands 
low density was found to be favorable for ferromagnetism.!! This tendency may be 
understood as a result of the diverging DOS at the zone boundary in one dimension. 
Since we have quite different results for d = 1 and d = oo for n < 1, a study of the 
ground states in two and three dimensions is desirable in order to answer the important 
question as to whether Hund coupling is effective in realizing ferromagnetism in low 
density systems like Ni. It should be noted that Ni has a fee lattice structure, and its 
DOS has a sharp peak near the edge, which is similar to the one-dimensional one. 



ELECTRONIC STATES IN THE DOUBLE EXCHANGE MODEL 

In this section we consider the so-called double exchange model (DEM), which is 
composed of electrons in a single conduction band and localized spins of magnitude S 
at all lattice sites. The electrons and localized spins interact through intraatomic Hund 
coupling. The DEM may be the simplest lattice model for electrons in manganites. If 
we assume a single orbital instead of the doubly degenerate e g orbitals and regard three 
electrons occupying t2 g orbitals as a localized spin in Mn 3+ ions, we obtain the DEM 
with S = 3/2. We may consider also the DEM for arbitrary S. For example, we may 
consider a level separation A between two orbitals in the doubly degenerate Hubbard 
model. If both U' and A are much greater than t, then the model reduces to the DEM 
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with 5 = 1/2 for n > 1. The DEM is described by the following Hamiltonian 

H = —t ( c L c j° + h - c -) ~ J J2 S *' w4t c jV, (4) 

(jj)Sn.n.,cr i,a,a' 

where s = ^ (<j x , a y , a z ) and cr Q denotes the Pauli matrix. The parameter J in (||) 
corresponds to 2 J in (|T|). Direct interactions between localized spins as well as the 
coupling between electronic and lattice degrees of freedom are neglected. For doped 
LaMnOs typicalvalues of the conduction band width 2W and J(25 + l)/2 are thought 
to be 1 ~ 2 eVE3 and 2 ~ 3 eV, respectively. For large JS the spin of an electron is 
always coupled parallel with the localized spin and forms a total spin of size 5 + 1/2. 
Since the original hopping term in @) conserves the spin of the electron, the hopping 
probability between two neighboring sites is effectively reduced when localized spins 
at these sites are not parallel to each other. The factor of the reduction is given by 
cos 9/2 if the localized spins are classical and make the angle 9 between them.0 In 
the paramagnetic state the localized spins are oriented randomly and as a result the 
conduction band is narrowed due to the reduction of the hopping integrals as well as 
decoherence effects due to scattering. Band narrowing increases the kinetic energy of 
the paramagnetic state and favors the ferromagnetic state. Ferromagnetism due to 
this "double exchange" mechanism was studied earlier O E3 The electronic states of the 
model with classical localized spins were studied by use of dynamical mean field theory 
by Furukawa.EJ Classical localized spins are not affected when they scatter conduction 
electrons. Quantum mechanical spins may be flipped during the scattering processes. 
The effect of this "spin exchange scattering" was studied earlier using the coherent 
potential approximation (CPA). The electronic states of the model with quantum 
spins in one and two dimensions were recently studied by using numerical methods 
extensively.^' 

The CPA theory mentioned above treated a single electron in a system with ran- 
domly oriented localized spins and did not take into account the presence of other 
electrons. As a result the theory is valid only in the low density limit, but in this limit 
it gives a qualitatively correct description of the change in the electronic states due to 
interactions with localized spins. The single conduction band is modified by the inter- 
actions and the density of states splits into two bands for JS J> W. The lower band 
corresponds to electronic states with electron spins parallel to the localized spins (we call 
them "parallel electrons") and the upper band to those of "antiparallel electrons". The 
relative weights of the lower and the upper bands are (5 + 1)/ (25 + 1) and 5/ (2S + 1), 
corresponding to the total spin of a site 5 + 1/2 and 5 — 1/2, respectively. If we naively 
consider these bands as a rigid one-particle density of states, the half-filled system (n = 
N e /N = 1) cannot be an insulator since the lower band is not full. Surely the half-filled 
system should be an insulator if JS 3> W. We need a theory which realizes the insu- 
lating half-filled system for JS ^> W in order to discuss the transport properties of the 
DEM. 

Let us first consider the atomic limit, i.e. the case with t = 0. The energy spectrum 
of the Green function in this limit is composed of four levels. The lowest level at 
u = — J(5 + l)/2 corresponds to the process of creating a parallel electron at a site 
which is already occupied by an antiparallel one. The second one at —JS/2 comes 
from creating a parallel electron at an empty site. The levels at J5/2 and J(5 + l)/2 
correspond to creating an antiparallel electron at an occupied and unoccupied site, 
respectively. The spectral weights of these four levels are (Sn — 2(Si ■ Sj))/(25 + 1), 
[(5 + l)(2-n) + 2(5 i -s,)]/(25 + l), [(5 + l)n + 2(S i -s i )]/(25 + l) and [5(2 -n) - 2(5; • 
Si)]/(2S + 1), respectively, where Sj denotes the electron spin operator at the i-th site. 
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When the hopping term is turned on, these discrete levels will broaden and compose 
four separate bands for t <C J. They broaden as t increases and will finally merge into a 
single band for f > J. We note that the first and third levels vanish in the low density 
limit (n = 0) and that the weights of the two bands of the spectrum in the earlier 
one-electron CPA theory reproduce those of the second and fourth levels correctly. As 
for the metal-insulator transition, we note that the two lower levels are fully occupied 
at n = 1 in the exact atomic limit Green function. Therefore, we can expect a Green 
function to reproduce the correct insulating behavior at n = 1 for JS W, if it 
reduces to the exact one in the atomic limit. We report below on the derivation of a 
Green function in the one-site approximation which gives the correct atomic limit and 
also reduces to the one-electron CPA result in the low density limit. We calculate the 
resistivity p in the framework of this approximation. A brief account of this work has 
already appeared and full details will be published elsewhere.^ 

In a single-site approximation the Green function G k {ui) is written as 



G k » = {G a (ujy L - (e k - J^))r, (5) 



where 



J (T {uj)=uj-Y I(T {lo)-G (T {lo)-\ (6) 

Here and Yi a {u)) denote the free band energy and the self energy, respectively. The 
local Green function G a (oj) is related to the self-energy through the DOS of the free 
band energy Do(efc) as 

J LO — h a {LU) — X 

In order to close the above equations we need another equation for G a (ui). We follow 
the equation of motion of the Green function along the line of Hubbard.C3 We can close 
the equations for arbitrary S in the paramagnetic state. In the paramagnetic state we 
can omit the spin suffices and the self-consistent equation for G(uj) is written as 

~ v (E(u) + aJ/2)(n/2) a + aJ/2(Sj • s_j) 

{ ' £^ (E{u) - aSJ/2) (E(u) + aJ(S + l)/2) ' 1 ' 

where E(oj) = uj — J(u>) and (n/2) a = 5 Q _ + an/2. The above set of equations reduces 
to that of the one-particle CPA at n = 0. We show in Fig. || the one-particle density 
of states for J = 4W and 5 = 3/2 obtained from above set of equations by using the 
elliptic DOS given by 

D Q (e) = 2/(ttW 2 )VW 2 - e 2 . (9) 

At n = the spectrum is composed of two bands centered at u = —JS/2 and J(5 + 
l)/2, respectively. As n increases the third band centered at uj = JS/2 emerges and 
correspondingly the weight of the band at lu = J(S + l)/2 decreases. There should be a 
fourth band at u = —J(S + l)/2 as well, but its weight is very small for JS ^ W, since 
{Si ■ Si) ~ nS/2 (In the calculation shown in Fig. |^ we approximated as (Si ■ Sj) = nS/2 
for simplicity). The third band grows with n and finally it takes over the second one at 
n = 1. At n = 1 the lower two bands (the lowest band in Fig. |5|) are completely filled 
and the system becomes a Mott insulator at T = as is expected. These three bands 
should be observable by photoemission experiments on manganites. We note that the 
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Figure 5: The density of states obtained from the approximate Green function in the 
paramagnetic state for J = AW and S = 3/2 at the filling n = 0, n = 0.25, n = 0.5, 
n = 0.75 and n = 1.0. 

position of the band at oj ~ JS/2 will be shifted to oj ~ JS'/2 + {/ in the presence of 
the intraatomic Coulomb repulsion U. For J = oo the lowest band has width 

2VF = 2WyJ(S + l-n/2)/(25 + l) (10) 

for the elliptic DOS. The band narrowing factor is a minimum at n = 1 with the value 
l/\/2 independent of <S. It should be noted that this factor may depend considerably 
on the choice of Dq as was found for one-particle CPA.H 

We calculate the DC resistivity p by using the above paramagnetic Green function in 
the Kubo formula. Vertex corrections do not enter in the calculation of the conductivity 
in the single-site approximation. The expression of the static resistivity at T = is 
simplified by assuming the cubic tight-binding form of the hopping term in H. We find 

P _1 = TIT / ede^e)D (e), (11) 
San J 

where d = Aj^(fj,) 2 \ e ^ =t and A^{p) is the spectral weight function of the Green func- 
tion Gfc (u>) at u) = fj,. This expression is evaluated using the elliptical approximation 
to the DOS, both in Dq(c) itself and in <p(e) via the Green function calculated above. In 
Fig. |6| we show the resistivity obtained for J = oo as a function of n for various values 
of S. We used 5A for the lattice constant a. Note that the correct insulating behavior 
is obtained for n = and n = 1. We find that p hardly depends on J for JS <; 5W. 
The resistivity p is of order of mficm for 0.1 ^ n ^ 0.9 and this is much smaller than 
typical experimental values for doped LaMn03 except for the case of Lai-^Sr^MnOs 
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Figure 6: Resistivity at T = (the paramagnetic state is assumed) is depicted versus 
n for J = oo. Data for S = 1/2 ~ 5/2 and oo obtained by using the elliptic DOS are 
shown. Resistivity increases with increasing S. 

with x(= 1 — n) ~ O.3MH The result shows also too weak a dependence on n compared 
to the experimental result. Furukawalll calculated p by using the Lorentzian DOS for 
S = 00 and obtained good agreement with experimental result at x = 0.2. We also 
calculated p by using the Lorentzian DOS defined as 

D (e) = W/[n(e 2 + W 2 )] (12) 

to calculate 0(e) but retaining the elliptical approximation to Dq(c) so that the integral 
in eq. (|ll]) converges. The results are shown Fig. [?]. They show much stronger depen- 
dence on n than those obtained for the elliptic DOS, and increase very rapidly when the 
fermi level approaches the band edge. The magnitude of p at x ~ 0.2 is several tens of 
niacin, which is of the same order as the experimental data. However, the elliptic DOS 
is considered to be more realistic than the Lorentzian one, since the Lorentzian DOS 
gives divergent second moments. Hence the good agreement with experiment obtained 
for the Lorentzian DOS may be an artifact. Our results suggest that scattering by 
random localized spins is not enough to explain the correct order of magnitude of the 
resistivity in manganites. The present result is obtained by treating dynamical aspects 
of the scatterings approximately. Effects of finite temperature and short-range corre- 
lations between localized spins are also neglected, since our calculation has assumed 
complete Fermi degeneracy and completely random configurations of localized spins. 
These effects may modify the above result to some extent, but we do not expect that 
they will change the order of magnitude of the resistivity. The ref ore some other effects 
should be taken into account to explain experimental results .0' e3 e3 

We also studied the magnetic properties of the system by using the Green function. 
The self-consistent equation for the general magnetization for 5 = 1/2 was obtained as 

<2 T M = (13) 
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Figure 7: Resistivity at T = for Lorentzian DOS. 

h± E^ a (u) (Ef(w)Ef(w) - J 2 /4) 

where E%(u) = uj—J C7 {uj)+aJ /A. We calculated the magnetic susceptibility x for J = oo 
by including a magnetic field and expanding (Sf + sf) about the paramagnetic state. 
We find that x never diverges at a finite temperature for any < n < 1, i.e. there is no 
ferromagnetic transition. At n = the correct Curie law \ = (9fJ>B) 2 S(S + 1)/(3/5bT) 
with 5 = 1/2 was obtained. On the other hand for n = 1 and J = oo x correctly obeys 
the Curie law with S = 1 at high temperatures but it obeys the law with S(S+1) = 2/15 
at low temperatures. 

In the DEM with J = oo the ground state is proven to be ferromagnetic in one 
dimension for any < n < 1. Ill At present no reliable study of the ground state phase 
diagram of the DEM with S = 1/2 seems to be available in higher dimensions than one. 
In three dimensions a high-temperature series expansion analysis suggests a finite Curie 
temperature for all electron density between and 1,H though the fully polarized state 
is not stable for 0.12 < n < CJ.45.e3 It is reasonable to expect that the ferromagnetic 
ground state is stable in three dimensions in some density region. Therefore we consider 
that our approximate Green function fails to reproduce the low temperature properties 
of the model correctly. It is known that the analogous CPA in the Hubbard model 
does not give ferromagnetism at any densitylil or the correct Curie law at n = 1.0 
Our present approximation suffers from a similar failure and an improved treatment is 
necessary in order to discuss the magnetic properties. We believe, however, that the 
present theory gives a qualitatively correct picture of the system in the paramagnetic 
state and that the paramagnetic resistivity obtained above is of the correct order of 
magnitude. We need an improved approximation to the Green function to study the 
magnetic properties of the system. A similar approach to that of KawabataEil for the 
strongly correlated Hubbard model might be useful. 
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SUMMARY 



We discussed two topics on the Hund coupling in lattice systems employing sim- 
plified models. First we examined the effectiveness of the Hund coupling in realizing 
ferromagnetism in the doubly degenerate Hubbard model. In quarter-filled systems 
the insulating ferromagnetic state accompanied by alternating orbital order was found 
stable. In more-than-quarter filling case metallic ferromagnetism is stabilized by the 
"double exchange mechanism". The above results are common to one and infinite di- 
mensions and we expect them to hold in general dimensions. In less-than-quarter filled 
systems the ferromagnetic ground state is stable in one dimension but not in infinite 
dimensions. To study this case in two and three dimensions is an interesting future 
problem. 

Secondly we examined the electronic states and the resistivity in the double exchange 
model by using the one-particle Green function. The splitting and narrowing of the one- 
particle spectrum due to the Hund coupling were clarified in the framework of a single- 
site approximation. The resistivity due to the scattering by random localized spins was 
shown to be too small to explain the experimental results of doped manganites. The 
present approximation failed to give the ferromagnetic state and we need an improved 
treatment to study the properties at low temperatures. 
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